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Abstract 

It has been recently pointed out that any primary galactic population 
of Weakly Interacting Massive Particles (WIMP) generates, through colli- 
sions with solar matter, a secondary population of "slow" WIMPs trapped 
in the inner solar system. We show that taking into account this "slow" 
solar-system population dramatically enhances the possibility to probe the 
existence of stable massive neutrinos (of a 4th generation) in underground 
neutrino experiments. Though neutrinos, with mass in the 45-90 GeV range, 
can only represent a sparse subdominant component of galactic cold dark 
matter, a combination of enhancement factors makes it possible to discrimi- 
nate their contribution to WIMP annihilation effects in the Earth. Our work 
suggests that a reanalysis of existing underground neutrino data should be 
able to bring extremely tight constraints on the possible existence of a stable 
massive 4th neutrino. 

1 Introduction 

The total number of quark-lepton families is not known for sure. Experi- 
mentally, three generations have been found at present. Theoretically, the 
number of generations could be larger. In superstring-inspired particle mod- 
els the number of generations is defined by some topological characteristics 
of the manifold of compactified additional dimensions, and any number of 
generations is a priori possible 

Another important prediction of almost all realistic superstring-inspired 
models is the existence of at least one additional U{1) gauge group in the 
low energy limit of the theory. Recently it was suggested to ascribe such 
a new U{1) gauge group to an additional, fourth, fermion generation only. In 
this case the new gauge group can remain unbroken. Such an unbroken f/(l) 
gauge group implies the existence of a strictly conserved charge, which, in 
turn, accounts for the stability of the lightest particle of the 4th generation. 
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and forbids any mixing with the other (usual) three generations. It will be 
assumed that the lightest particle of the 4th generation is its neutrino. 

The direct search for new generation fermions on accelerators leads to a 
lower bound on their masses in the range 50-100 GeV. The mass of a 4th 
generation neutrino is restricted by the measurement of the width of the 
Z-boson to nil, > rnz/2 ^ 45 GeV if it is a (quasi) stable Dirac neutrino (for 
a Majorana one the restriction is slightly lower; for an unstable one the lower 
limit is about 90 GeV). Another possibility to search for a fourth generation 
neutrino, with mass mj, > mz/2 ^ 45 GeV, at accelerators was suggested in 
||2|, 0]. A detailed analysis of the data on the parameters of the Standard 
Model, accounting for the possible contributions of virtual new generation 
fermions, allows for the existence of an additional generation if the mass of 
the new (fourth) neutrino is about 50 GeV Q and if the masses of the other 
fermions of the new generation exceed 100 GeV. 

The existence of new generation fermions in the Universe can lead to 
many observable astrophysical effects. This makes the appropriate cosmo- 
logical and astrophysical analysis an important tool for probing the possible 
existence and properties of such 4th generation fermions. In particular, a 
new neutrino, being plausibly the lightest of its generation, and thus possi- 
bly stable, is of the most interest in such an analysis. 

As was found long ago [|], the existence of a heavy Dirac neutrino is com- 
patible with the upper limit on the total density of the Universe if its mass 
exceeds 2 GeV. Indeed, for a neutrino mass greater than 3 MeV, neutrino 
annihilation through weak interactions reduces their cosmological concentra- 
tion at freeze out. A larger neutrino mass corresponds to a larger annihi- 
lation cross section, and therefore to a smaller relic heavy neutrino density 

= Pu/ Perineal (scc Fig. 2 bclow). Morcovcr, if m,^ = mz/2 the Z-boson 
annihilation resonance leads to a dip in the relic density. For instance, if 
mu = 50 GeV, close to the Z-boson resonance dip, fl„ ~ lO"'^. Such a rare 
population of neutrinos does not play any significant dynamical cosmological 
role as a Cold Dark Matter (CDM) contributor. Nevertheless a more refined 
astrophysical analysis 0, Q showed that the effects of such rare, massive 
neutrinos can be accessible to some experimental searches and/or astrophys- 
ical observations. In the present paper, we shall mainly consider the mass 
range 45 < rriu < 90 GeV, corresponding to such a sparse population of heavy 
neutrinos. 

Heavy neutrinos, as any form of CDM, must be concentrated in galaxies. 
The ratio of galactic neutrino density to the mean cosmological density is a 
model-dependent parameter (denoted here as ^) which is strongly sensitive 
to the details of galactic halo formation. In this work we shall assume that 
the concentration factor ^ is the same for massive neutrinos and for the 
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dominant contributor to CDM. It is usual to estimate ^ by taking the ratio 
between a local (i.e, in the vicinity of the Solar system) density of CDM 
equal to 0.3GeV/cm~^, and a mean cosmological density corresponding to 
^CDM = 0.3. This leads to the "standard" estimate: ^ = 2 x 10^. In most of 
this paper, we shall assume this standard value (except, when explicitly said 
otherwize). 

In previous papers m, [|] we already analyzed several astrophysical con- 
sequences of such a sparse galactic population of heavy neutrinos. In partic- 
ular we studied [0] the (ordinary) neutrino fluxes emitted by the annihilation 
of heavy neutrino-antineutrino pairs accumulated in the core of the Earth. 
[Such an accumulation takes place for many types of weakly interacting mas- 
sive particles (WIMP) ||.] 

Recently, it was pointed out Q that this annihilation flux can be strongly 
enhanced by the existence of a "slow" Solar-system population of WIMP's, 
trapped in the gravitational held of the Solar system by an initial inelastic 
interaction with the Sun |[lO|| . The aim of the present paper is to analyze in 
detail, in the case of a heavy neutrino WIMP, the density of this new "slow" 
Solar-system population, and its enhancement effect on the annihilation flux 
from the core of the Earth. We shall show that the existence of the slow 
population qualitatively improves the sensitivity of underground neutrino 
data to the effects of 4th neutrino annihilation, and makes these data a 
signiflcant probe of the existence of a 4th neutrino. 

2 Number density of the "slow" heavy neutrino 
population 

Let us recall (from ||l^) that the "slow" population of WIMPs is generated 
by inelastic collisions of incident primary ("fast") galactic WIMPs with nuclei 
in the outer layers of the Sun. A fraction of these WIMPs are scattered, 
by the collision, on orbits that "graze" the surface of the Sun, and which 
evolve, under the subsequent perturbing gravitational influence of planets, 
onto orbits that do not penetrate the Sun. This allows these WIMPs to 
survive for a long time in the Solar system. This population has a much 
lower typical velocity (~30 km/sec) than the incident galactic one. This fact 
amplifles the probability of their capture by the Earth (we will further refer 
to this population as "the slow component"). The number density (in the 
vicinity of the Earth) of this slow component of heavy neutrinos, risiow, has 
been derived in Ref.||lO| and reads: 
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Here /(f oo) is the angular average of the velocity distribution of galactic 
WIMPs. [Our notation differs slightly from Ref. ||10| in that we factor 
out the number density Ugai from the phase-space distribution function, and 
delete the bar over /(foo) indicating the angular average.] The parameter 
Vq is the characteristic velocity entering the (assumed) Maxwellian velocity 
distribution /(woo). Throughout the paper the index oo refers to quantities 
far from the Sun or (depending on the context) the Earth. The meaning of 
the other quantities entering the equations above is: ttia is the mass of the 
nucleus A, denotes its mass fraction in the Sun, cta is the cross-section of 
neutrino-nucleus scattering in the point-like approximation. The form-factor 
F{voo) takes into account the effect linked to the extension of the nucleus. 
More about cross-section and form-factor below. The quantity gtot (which 
has dimensions [crosssection] x [mass]~^ should be expressed in GeV'^. Here 
and below we use the notation 

nil, denoting the neutrino mass. The upper limit, fmax? of the integral 



above comes from kinematics and is read from Eq.(2.16) of [10 



4m^mAV^sc{Rs) - {my + tuaYcx 



\my - mA\ 

in which we used the following numerical estimates: 



(3) 



vlXRs = 0.907i?5) = (644km/s)2, (4) 

where Rs denotes the Solar radius, and a = GMs/a ~ (30km/s)^, a 
denoting the semi-major axis of the WIMP orbit. As said above, we take 
a WIMP velocity distribution (far from the Sun) which is Maxwellian, i.e. 
(after angular averaging) 
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Here vs denotes the velocity of the Sun relative to the gas of galactic 
neutrinos. In the estimates below, we took vs = 200 km/s. We also estimated 
the effect of varying vs between zero and 200 km/s and found that it changed 
only insignificantly the results given below. The parameter vq denotes the 
characteristic velocity dispersion of the Maxwell distribution (i.e. it measures 
the "temperature" of this distribution). In our estimates below we shall use 
the value vq = 220 km/s. For comparison, we shall also mention the results 
for fo = 266 km/s, which corresponds to a mean velocity of 300 km/s. 

The data about the chemical composition of Solar matter was taken as 
in Ref. ||T^ (i.e., a combination of the values given in and ||12|)- We also 
follow Ref. [|l^] in taking nuclei form-factor of the form: 



exp 



Q{v) 

Qa 



(6) 



Q 



where Q 



2 is the transferred energy from neutrino to nucleus, 
: 10-13 cm (0.3 + 0.91( 



^ GeV ) 



The total u — A cross-section for a point-like nucleus a a is approximated 
by a purely spin independent coupling, except in the case of hydrogen. Since 
most contributing nuclei have a first excited energetic level higher than the 
characteristic energy tranfer Q the collision is assumed to be elastic. This 
leads to a cross-section of the form 



2 2 

mi,m\ 



'^A = ^ , T"^2 (^v + 3FX), (7) 
Stt [rriu + mA) ^ ' 

where Gp is the Fermi constant, Fy = —{A — 2(1 — 2siia^ 9w)Z) is the 
vector form- factor (all g^-dependence is taken out into F{v)), 9w is the Wein- 
berg angle, and Fa is the axial form-factor, which is taken into account only 
in the case of hydrogen (hence the Kronecker symbol Sa)- 

After substituting all these quantities in Eq. (1) we obtain the dependence 
on rrii, of the density enhancement ratio 6e = ^^^^^ showed in Fig.l. In Table 
1, we fix rriu = 50 GeV, and study the sensitivity of the density enhancement 
Se to the choices of the two velocity parameters entering our study: Vq and 
Vs- 





vo = 220 km/s 


Vq = 266 km/s 


vs = 


0.728 


0.529 


Vs = 200 km/s 


0.475 


0.388 



Table 1: Density enhancement 6e 
ters. 



nsiow/ngai for different velocity parame- 
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Figure 1: Density enhancement 6e = nsiow/ngai as a function of the neutrino 
mass nil,. 



The crucial new result to note from Fig. 1 and Table 1 is that, contrary 
to the case of a generic neutralino WIMP where the density enhancement 
(near the Earth) due to the slow population, 6e = nsiow/ngai, was typically 
of the order of a few percent, here this density enhancement is of the order of 
50 %. This shows that it is crucial to take into account the existence of the 
slow population to estimate the detectability of a 4th neutrino . In addition 
to this significant enhancement in the local density of WIMPs, we turn now 
to the estimate of the capture rate of the WIMPs by the Earth which is 
further amplified by the fact that slow WIMPs (with velocity ~ 30 km/s) 
are much more easily captured by the Earth gravitational field than galactic 
ones (with velocity ~ 220 km/s). 

3 Capture rate of massive neutrinos by the Earth 

The neutrinos interacting with the Earth matter can be captured by the 
Earth gravitational potential well if they lose enough energy in the colli- 
sion with a nucleus in the Earth. The "slow" neutrinos, having an inci- 
dent velocity already comparable with the escape velocity from the Earth ( 
f g ~ 10 km/ s) , have a much greater probability to be captured than galactic 
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ones. It is therefore quite important to take their effect into account. 

Among all nuclei present in the Earth the iron nuclei play the main role 
in the capture process for the neutrino mass range ( 45 GeV < my < 90 GeV) 
that interests us. 

The capture rate of WIMPs by the Earth has been studied in detail by 
Gould [|l^] (see also Q). In the present study, we estimated the capture 
rate by starting from Eq.(5.7) and Eq.(5.8) of 0, and by approximating the 
result in the following way. We replaced the integrals over the volume of the 
Earth of the type / d^xnA{r)f{r) by NA{f{r)), where Na is the total number 
of nuclei of type A, and where {f{r)) is a suitable (approximate) average of 
/(r). Actually, we consider only the capture by iron. As iron is concentrated 
in the core of the Earth we can use the average of the squared escape velocity 
given in [|^: {v^^^) = 1.6 x (11.2 km/s)^. We assume that the fraction of the 



mass of the Earth in iron is 20 %. Finally, our simplified capture rate reads: 



n ^Qfs f 
P-rrii, Jo 







exp ^ 





(3-171^ 



exp 



2Q 



Fe 



(8) 



The neutrino distribution function driu to be inserted in Eq.(8) depends 
on whether we consider galactic or slow neutrinos. For galactic neutrinos it 
is dngai = Ugai f{v)d^v. For slow neutrinos we followed Eq.(3.15) of (with 
A = 1 ; £ = 0.18377), i.e, explicitly. 



X 



(^l^slow slow J slow \^ slow slow f^slow^^ /Q 2 2^06 ^ 

[6 X — V ) ■ 

(^^Jv^ -0.816 X vie {v - 0.904 x ve) - \/v^ - 1.184 x vl9{v - 1.089 x ve)^ dv , (9) 



where C = 91.775 ^ is a normalization factor, 9{x) is the step function, 
Ve = 30km/s is the velocity of the Earth orbital motion, and where 

^slow 1^ 

given by Eq.(l) above. 

Another crucial parameter for our problem is the value of the incident 
galactic neutrino number density itself. Contrary to the case usually consid- 
ered for WIMP capture, we cannot assume here that the sparse 4th neutrino 
population is the principal constituent of Cold Dark Matter, i.e we cannot as- 
sume that its galactic density has the "standard" value Pcdm ~ 0.3 GeV/cm^. 
Instead, we need to estimate what is the relic density of 4th neutrinos left 
over from the Big Bang, and to multiply it by an estimate of the typical 
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density enhancement, in a galaxy. We can then write the galactic neutrino 
number density (as a function of the neutrino mass) as: 



n(mi,) = irireiirriu). (10) 

Here ^ denotes the parameter of neutrino clustering in the Galaxy, i.e. 
the ratio of the local neutrino density to the relic one, rireiijny). As already 
mentioned above, in our estimates we take the standard value ^ = 2 x 10^. 

As for the relic number density (of neutrinos and antineutrinos) , we take 
the estimate derived in [Q, namely: 

, , 4.2x10"^^ , ( gs , n^\ . frn^\ 

nrei[mu) = —= -, ^ 40 + In ——mpm^ {(TannP} - In J— ]xn^ 



fglmpniy {aannP) \ Vv^ / V 

(11) 

Here ~ 80 is the number of effective degrees of freedom at the mo- 
ment of neutrino freeze-out, (yf^ = 2 is the number of neutrino spin states, 

n. 
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411 cm~^ is the present number density of relic photons, {(JannP) is 
the thermally averaged annihilation cross-section (multiplied by the dimen- 
sionless neutrino relative velocity) at freeze out. The value of the last log- 
arithm entering the above result can be taken to be 1.6 in the mass range 
45 GeV < m,^ < 90 GeV that we shall consider. In this mass range, the 
dominant annihilation channel for neutrinos is the channel involving one in- 
termediate Z-boson. In this case the value {(JannP) may be taken (for the 
non-relativistic case) in the form 

—4 2 

(aannP) = a.„n/3*V2 = Br'^Z ^e^e)|^ (4^2_^|)2 + ^|r| ^' 

(12) 

where Br{Z — > Uei'e) ~ 0.0667 is the branching ratio of the decay of the 
Z-boson into two electron neutrinos, mz and Tz are the mass and width 



of the Z-boson, and g = ^aV^Gf rriz is the dimensionless constant of weak 
interaction. Note that, if we were to consider larger neutrino masses (m^ > 90 
GeV), the channel w W^W~ would start to dominate and the cross- 
section would increase with rrij^. 

The relic neutrino density in units of the critical density is shown in Fig. 2. 
Note in particular that, when niy = 50 GeV, = 1.7 x 10~^ (which is indeed 
much smaller than the usually considered WIMP relic densities). Still in this 
numerical example of mj, = 50 GeV, by multiplying this relic density by the 
amplification factor ^ = 2 x 10^ mentioned above, we find a local, galactic 
number density of massive neutrinos equal to about 0.2MeV/cm^. This is 
indeed much smaller than the "standard" value for usual WIMPs. 
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Figure 2: The density of primordial neutrinos in units of the critical density. 

After their capture, neutrinos (and antineutrinos) are thermalized, set- 
tling in the Earth core, and neutrino-antineutrino pairs start to annihilate. If 
the equilibrium between capture and annihilation is reached, the annihilation 
rate is equal to the capture rate. Let us estimate the capture rate (Req) for 
which the equilibrium between capture and annihilation inside the Earth is 
reached during the lifetime of the Earth. The variation of the number of 
accumulated neutrinos and antineutrinos satisfies 



where Reap and Ram denote the capture and annihilation rates respec- 
tively. We have 



Here the factor 2 accounts for the fact that in one annihilation a pair 
of neutrinos disappears; the factor 1/4 comes from the fact that n denotes 
the total (neutrinos plus antineutrinos) neutrino density within the ther- 
malized core of the Earth; in {aannVrei), Vrei denotes the relative velocity of 
thermalized neutrinos. A rough estimate of Rann is obtained by assuming a 
homogeneous distribution of thermalized neutrinos within a certain volume 
Vtherm- Wc dcfiuc Vtherm as the volumc bouudcd by the radius which can 
be reached by a particle of kinetic energy Ttherm in the center of the Earth 
freely moving in the Earth potential well. The typical neutrino energy Ttherm 



dN 
~dt 



'Cap 



(13) 




(14) 
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is supposed to correspond to the temperature of the Earth core. The lat- 
ter is not known exactly, but is around 10000 K. We can therefore consider 
Ttherm = 1 cV as fiducial value. This yields: 

^therm — o'^^Rtherm^ Rtherm — \ J!^^ ^ ^ Re^ (15) 

-J V J-escO PEcore 

2 

where Tesco = "^'^"^'"'^ , (pe) = 5.5g/cm^ and pEcore ~ 12g/cm2 are the 
mean and core Earth densities, and Re is the Earth radius. In this approxi- 
mation we have 

Rann = 7^ {(^annVrel) =CN . (16) 

^ therm 

It is then easily checked that, neglecting the (slow, linear) variation with 
time of Reap, the solution for the time variation of N, or equivalently Rann, 
is Rann — -Reap tanh^ (^yJCRcapt^- The "equilibrium" rate Req, for which a 
balance between capture and annihilation within the Earth lifetime is 
maintained, is therefore estimated to be Req = l/{Ct'^). 

Substituting all the factors, we finally obtain 



Here the mass of the neutrino, as well as the mass and the width of the 
Z-boson are in GeV. For numerical estimates, we took = 2, and (as 

said above) Ttherm = 1 eV. 

Having in hands, ready for comparison, the "equilibrium" rate Req, we 
now come back to the estimate of the capture rate, as given by Eqs.(8)-(11) 
above. Note that the capture rate of both galactic and slow neutrinos are 
proportional to the assumed value of the concentration factor ^. 

We plot in Fig. 3 the total capture rate of neutrinos (summed over slow 
and galactic ones), as a function of the neutrino mass, for ^ = 2 x 10^. We 
also indicate what would be the capture rate if one included only galactic 
neutrinos. 

The dip at rrij, ^ 45 GeV reflects the resonant annihilation dip in primor- 
dial neutrino density (see Fig. 2). The peak in the galactic neutrino capture 
rate at rrii, ~ 55 GeV comes from the well known fact that when the WIMP 
mass is equal to the mass of iron nuclei their collision is efficient in slow- 
ing down the WIMPs. [This peak is shifted from the iron mass niFe — 52 
GeV because of the rising relic neutrino density factor.] In the case of the 
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Figure 3: Total capture rate of slow and galactic neutrinos (solid line), 
compared to the capture rate of galactic neutrinos only (dashed line). 
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Figure 4: Ratio of the slow neutrino capture rate to galactic one. 
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slow WIMPs (which lose their kinetic energy during collisions much more 
effciently than the fast galactic WIMPs) this "iron resonance" effect is spread 
over a large range of masses and does not show up as a peak. 

In Fig. 4 we present the ratio of slow neutrino capture rate to the galactic 
one. This ratio, obviously, is independent of galactic neutrino density. The 
dip around rrii, = mpe = 52 GeV in Fig. 4 is due to the "iron resonance" peak 
of the galactic neutrino capture. 

An important conclusion from Fig. 3 is that, in the considered mass range, 
the equilibrium between capture (Reap) and annihilation rates (Rann) inside 
the Earth is established on a time scale smaller than the Earth lifetime. 
Therefore the annihilation rate can be simply taken as being equal to the 
(above computed) capture rate: Ram = Reap- 
To complete the results indicated in Fig. 3, we give in Table 2 the galactic 
and slow neutrino capture rates (in units of 10^^ s~^) for different velocity 
parameters, and for the special value rrii, = 50 GeV (and for the standard 
value ^ = 2 X 10^). Note that the equilibrium capture rate for nii, = 50 GeV 
amounts to 0.53 x 10^^ s"^ 





vo = 220km/s 


Vq = 266 km/s 


galactic 


slow 


galactic 


slow 


vs = 


3.86 


18.6 


2.88 


13.50 


vs = 200km/s 


2.62 


12.1 


2.14 


9.88 



Table 2: Capture rates in units of 10 s of galactic and slow neutrinos 
when nil, = 50 GeV, for different velocity parameters. 

One can see that our results are robust under varying the velocity param- 
eters (retaining the Maxwellian form of the velocity distribution). 

In the estimates above, we have not taken into account the effect of the so- 
called "Sakharov enhancement" of the annihilation rate due to the Coulomb- 
like attraction mediated by the new U{1) interaction ascribed to neutrinos 
and antineutrinos of the 4th generation [|l^]. This interaction has two effects: 
on one hand it reduces the relic neutrino density (by about 10% if ay = 
1/137, and 25% if ay = 1/50), but on the other hand it strongly increases 
the annihilation rate in the Earth. The overall effect is to strongly reduce 
the equilibration time scale between capture and annihilation in the Earth. 
Finally, if we were to take into account the "Sakharov enhancement" we 
would be in the conditions of equilibrium for the whole neutrino mass range 
considered and for all the acceptable magnitudes of parameter ^. 
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4 Conclusions 



In the present paper we studied the capture by the Earth, and the anni- 
hilation in the Earth core, of hypothetical fourth neutrinos (candidate to 
a sparse sub-dominant component of galactic CDM). We took into account 
not only the primary "fast" population of neutrinos, but also the recently 



pointed out secondary "slow" population ||T0|. It was found that the account 



of the slow component is crucially important in the considered neutrino mass 
range, 45 < rrii, < 90 GeV. Indeed, the contribution to capture of the slow 
component is larger by up to the two orders of magnitude than the one of 
the galactic component (Fig. 4). 

These results suggest the crucial significance of underground experiments 
(AMANDA, Super-Kamiokande, Baksan and others) for testing the fourth 
neutrino hypothesis. For example, Ref. ||14| has derived the constraint on 
WIMP annihilation in the Earth, from AMANDA data, under some assump- 
tions on the annihilation channels, and for a WIMP mass larger than 100 
GeV. Making a rough extrapolation of this constraint down to a mass of 
about 50 GeV one finds a potential sensitivity of the already existing under- 
ground neutrino data to the annihilation of a 4th neutrino in the Earth, for 
almost all the considered interval of neutrino mass. Note, that the presence 
of the slow component of a 4th neutrino plays a crucial role in this poten- 
tial sensitivity. Of course, this example can serve only as an illustration, 
since special analysis of the data in the framework of the hypothesis of a 4th 
neutrino is needed. 

The possibility to distinguish, in underground neutrino experiments, the 
contribution to annihilation effects of a sparse component of 4th neutrino 
from the contribution of other WIMPs (presumably dominating in the galac- 
tic CDM), results from the combination of several factors. The neutrino cap- 
ture in the Earth is facilitated by the relatively large weak interaction cross 
section of a massive neutrino and by the kinematic enhancement of neutrino 
momentum losses in collisions with iron nuclei. The neutrino annihilation 
effects in the Earth are strengthened by the relatively large neutrino weak 
annihilation cross section near Z-boson resonance (which is further strongly 
enhanced by the Coulomb like effect of the new long range interaction), 
and by the existence of a monochromatic neutrino-antineutrino annihilation 
channel, specific to a 4th neutrino. 

The presence of a slow component qualitatively enhances these factors. 
The slow component increases by up to 50% the number density of 4th neu- 
trinos near the Earth. Owing to their order-of-magnitude smaller mean ve- 
locity, the slow neutrinos are more effectively captured by the Earth (by up 
to two orders of magnitude) than the galactic ones. In the slow component 
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capture the kinematic peak of iron nuclei capture is spread over the whole 
considered neutrino mass interval, making it accessible to experimental test. 
The establishment of kinetic equilibrium between neutrino capture and an- 
nihilation in Earth makes the predicted annihilation fluxes insensitive to the 
details of captured neutrino distribution. As a result, the account of the slow 
neutrino component makes the hypothesis of a stable massive 4th neutrino 
accessible to underground neutrino experimental tests even under the most 
unfavourable astrophysical conditions. 

To conclude, our work shows that it is important to analyze existing (and 
future) underground neutrino data with the view of probing the hypothet- 
ical existence of a stable fourth generation neutrino with a mass about 50 
GeV. The analysis of the data of MACRO, AMANDA, Kamiokande, Bak- 
san and/or Super-Kamiokande is expected to provide an important probe 
of (and probably stringent constraints on) this hypothesis, especially in the 
case where one wishes to explain the DAMA event rate by assuming heavy 
neutrinos. The existence of a slow 4th neutrino component is crucial in such 
an analysis because it strongly enhances the underground neutrino flux ex- 
pected from 4th neutrino annihilation in Earth. This analysis can be viewed 
as a modest step towards the study of heterotic string phenomenology which 
generically leads to the prediction of an additional U{1), which, in turn, 
provides a motivation for considering a stable 4th generation neutrino. 
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